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ABSTRACT 

While both isothermal and NFW-based mass models for galaxy clusters are widely adopted in 
strong- lensing studies, they cannot easily be distinguished based solely on observed positions of arcs 
and arclets. We compare the magnifications predicted for giant arcs obtained from isothermal and 
NFW profiles, taking axially-symmetric and asymmetric mass distributions into account. We find 
that arc magnifications can differ strongly between the two types of density profiles even if the image 
morphology is well reproduced. Magnifications by lenses with NFW density profiles are usually larger 
than those for lenses with singular or nearly singular isothermal density profiles, unless the latter have 
large cores. Asymmetries play an important role. We illustrate our results with the two well-studied 
clusters MS 2137 and A 370. We confirm earlier results showing that both isothermal and NFW mass 
models can very well reproduce the observed arcs, radial arcs and other arclets. While the mass model 
for MS 2137 is not very well constrained, the two types of mass models produce strongly differing 
critical curves and caustics for A 370. We find that the NFW mass profile is preferred for A 370. We 
identify new arclet candidates in the field of A 370. Redshift estimates allowed by the lens model 
are consistent with results in the literature, where available. Three newly found counter-images are 
suggested to arise from an active, dust-enshrouded star- forming galaxy at z « 1.1. 
Subject headings: dark matter - gravitational lensing - galaxies: clusters: individual (MS 2137, A 370) 



1. INTRODUCTION 

Strong lensing by galaxy clusters, giving rise to 
strongly distorted arc-like images, plays an important 
role in studies of the mass distributions of galaxy clus- 
ters. It is highly non-linear, weakly affected by baryonic 
physics, and the data reduction is comparatively simple 
because it focuses on the positions of arcs and arclets 
(e.g. Leonard et al 2007; Limousin et al 2006; Halkola, 
Seitz & Pannella 2006; Zekser et al 2006; Covone et 
al 2006; Comerford et al 2006; Broadhurst et al 2005; 
Gavazzi 2005; Gavazzi et al 2003; Smith et al 2005; Kneib 
et al 2003; Kneib et al 1993; Mellier, Fort & Kneib 1993; 
Hammer 1991). Cosmological applications of strong lens- 
ing range from the use of clusters as gravitational tele- 
scopes (Frye et al 2007; Ofek et al 2006; Chary, Stern & 
Eisenhardt 2005; Kneib et al 2004; Ellis et al 2001) to 
cosmological constraints based on arc statistics (Fedeli 
et al 2008; Hamana et al 2005; Li et al 2005; Meneghetti 
et al. 2005a,b; Puchwein et al 2005; Dalai & Holder, 
2004; Wambsganss et al 2004; Bartelmann et al 1998). 
Currently, about 100 clusters are known to contain large 
arcs, including those newly discovered by Bolton et al 
(2006) in the Sloan Digital Sky Survey data. 
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One way of investigating cluster mass distributions 
with strong lensing starts from an assumed mass profile, 
which is commonly taken to be the well-known singular 
or non-singular isothermal sphere (SIS or NIS, respec- 
tively) 

1 1 

Psis oc — or pnis oc 2 , (1) 

T T ~\- T Q 

where r is the radius from the cluster centre and r c the 
core radius. 

Another mass model commonly used for strong lensing 
is based on the high-resolution numerical simulations of 
dark-matter halos in the ACDM framework by Navarro, 
Frenk & White (1996, 1997; hereafter NFW), who found 
the density profile 
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parameterized by the scale radius r s . While simulations 
and X-ray observations of galaxy clusters both support 
the NFW profile and applications to strong-lensing stud- 
ies have been successful (Zekser et al 2006; Comerford 
et al 2006; Kneib et al 2003), the SIS and NIS mod- 
els can also reproduce observed arcs remarkably well in 
many clusters (e.g. Kneib et al 1993; Mellier et al 1993; 
Miradal-Escude 1995; Gavazzi et al 2003; Gavazzi 2005). 

So far, conclusive evidence that the density profiles of 
real and simulated galaxy clusters agree with each other 
has not been found. Unfortunately, it is virtually impos- 
sible to distinguish isothermal and NFW density profiles 
based solely on the positions of observed arcs because 
of the near-degeneracy between these two density pro- 
files close to the scale radius. If radial arcs or multi- 
ple tangential arcs at different distances from the cluster 
centre are available, the situation improves substantially 
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because they probe the density profile at different radii 
(Gavazzi et al 2003; Gavazzi 2005; Tu et al 2008). While 
the positions of arcs and arclets are determined by the 
first derivative of the lensing potential, their magnifica- 
tions and distortions depend on its second derivatives. 
Different models reproducing the image geometry can 
thus yield substantially different image magnifications. 

The scientific goal of this paper is to investigate how 
relative arc and arclet magnifications can be used to 
break degeneracies between different mass models for 
strongly lensing clusters when combined with positional 
and photometric information. 

Two well-studied clusters, A 370 and MS 2137, are 
used for illustration. The code for fitting lens models 
is based on Comerford et al (2006) and will be briefly 
described below. Wherever needed, we adopt the "con- 
cordance" flat cosmology with fl m = 0.3, J7a = 0.7 and 
Hubble constant h — 0.7. 

The paper is structured as follows. Starting from a 
one-component lens model, we compare in Sect. 2 the 
predicted magnifications of the resulting arcs for the SIS, 
NIS, and NFW profiles, respectively, starting from axial 
symmetry and taking deviations therefrom into account. 
The results are easily extrapolated to lens systems with 
multiple components since the distance between two lens 
components is typically larger than their core or scale 
radii, r c or r s . In Sect. 3, earlier work on MS 2137 
and A 370 is briefly summarized. We adapt different 
mass models to these clusters, derive the arc magnifi- 
cations and compare with observations. Based on the 
high-quality images observed, we study the arclets in the 
field of A 370 in detail. We present our conclusions in 
Sect. 4. 

2. THEORETICAL PERSPECTIVES 

2.1. Basic lensing notation 

We first outline the notation to be used below. As 
usual, lenses are assumed to be sheet-like perpendicular 
to the line-of-sight. Then, the lens mapping from the 
lens to the source plane can be written as 

S = L - W , (3) 

and the local distortion of the lensed images is deter- 
mined by the Jacobian matrix 



where (Si, $2) and (Li,^) are dimensionless coordi- 
nates in the source and lens planes, respectively, and ip 
is the effective lensing potential, which characterizes the 
lensing mapping completely for any given lens model. 
The magnification at a given position L is 

^dibr (5) 

Giant arcs appear close to the critical curves when their 
sources are close enough to the corresponding caustics. 
Since giant arcs are the most pronounced structures due 
to highly non-linear strong lensing, they provide sensitive 
constraints on the lens models. 

For simplicity, we adopt in the following subsections 
a one-component lens model to illustrate the predicted 
magnifications of expected giant arcs for the SIS, NIS, 
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Fig. 1. — The predicted gradient (ddet^/ds)|t of the Jacobian 
determinant at the position of the tangential critical curve is shown 
as a function of the dimensionless radius s with the solid, dotted 
and dashed lines for the SIS, NIS, and NFW density profiles, re- 
spectively. Dotted lines from top to bottom show the results for 
the NIS model with r B /r c = 50, 10, 5, respectively, compared to the 
scale radius of the NFW model (see the text for details). 

and NFW profiles, respectively. Although real lenses typ- 
ically have multiple components, the main conclusions on 
the magnification will remain valid because they are due 
to generic lensing behaviour near critical curves or caus- 
tics. 

2.2. Axial symmetry 

If the lens is axially symmetric, i.e. the ellipticity e = 
0, we can obtain the radial gradient (ddet*4/ds)| t of 
the lens mapping Eq. (g]) at the radius r t where giant 
tangential arcs appear. There, det „4 = 0. We introduce 

s = r/r s (6) 

throughout this section for easier comparison. We plot 
the gradient (d det _A/ds)| t as a function of s in Fig. [TJ 
where the solid, dotted and dashed lines show the results 
for SIS, NIS, and NFW profiles, respectively. 

The radial gradient of the Jacobian determinant at the 
tangential critical curve is smaller for NIS profiles with 
different core radii than for the SIS profile, since the NIS 
density profile is always flatter than that of the SIS, thus 
the projected mass distribution is flatter for the NIS. 
With increasing core radii in NIS profiles, i.e. decreasing 
r s /r c , the gradient decreases due to the flattening of the 
density profiles. If the giant arc is located at a radius 
r t ^ r s , i-e., s > 1, the gradient is larger for a NFW 
than for the SIS and NIS profiles because it is steepens 
towards oc r~ 3 outside the scale radius. If the giant arc 
appears in the inner region, i.e. r t < r s , the gradient for a 
SIS or a NIS lens with a small core (r s /r c > 10) is larger 
than that for a NFW lens, while the gradient for a NIS 
profile with a relative large core (r s /r c < 10) is smaller 
than for a NFW profile due to its flatter density profile. 

It is interesting to study the corresponding magnifica- 
tions of the giant arcs produced by these different mass 
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profiles. The determinant det A of the lens mapping in 
Eq. ((4|) at a distance 5r from the critical curve can be 
written as 



det„4 = detA| t + 



ddetA 



dr 



5r = 



ddet.A 



dr 



(7) 



where the subscript t again denotes the position of the 
tangential critical curve where det„4 = 0. If the width of 
the arc is W , its total magnification p, t can be expressed 

by 



J dSj = J dSj 
J dS s J dSi/fii 
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If an observed giant arc appears at a radius r t ~ r s ~ 15" 
from the cluster centre with a width of W ~ 3", say, its 
estimated magnification can reach ~ 100. 

The predicted total magnification p, t scaled by (rt/W) 2 
for giant arcs as a function of their distance s from the 
lens centre is plotted in Fig. [2] with the same meaning 
of the line types as in Fig. [1] The figure shows that the 
expected magnification of an arc is always lower for a SIS 
than for a NIS model, given the observed positions of gi- 
ant arcs and their widths. This is because the SIS profile 
is always steeper than the NIS profile. If the arc appears 
a t r t > r s , the predicted magnification for a NFW profile 
is the lowest due to the steepest density distribution. If 
the arc appears at r t < r s , the expected magnification for 
the SIS or the NIS models with a small core (r s /r c > 10) 
is lower than that for a NFW model. Only when the 
core of the NIS profile is large (r s /r c < 10), compara- 
ble to the scale radius of a NFW profile, the expected 
magnification is larger than that for a NFW model since 
its density profile will be significantly flatter in the inner 
region. The conclusions can thus be well understood in 
terms of the flatness of the density profile near the arc 
locations. See also the discussion in Meneghetti et al 
(2005b). 

2.3. Asymmetry 

Since the gravitational potentials of real lensing clus- 
ters must be asymmetric, as was already pointed out 
by Grossman & Narayan (1988), Kovner (1989), Mellier 
et al (1993) and Miralda-Escude (1995), we now turn 
to asymmetric models, which play an important role in 
cluster- lensing studies (Bartelmann & Meneghetti 2004) . 
As in Bartelmann & Meneghetti (2004), we define the 
surface-mass distribution of the cluster in terms of the 
projected elliptical radius 



r c = [(rcos6») 2 (l - e) + (rsin6») 2 /(l - e)] 



1/2 



(9) 



Here, 9 is the position angle counterclockwise from the 
+y axis. The ellipticity 

e = ^ (10) 
a + b 

of an ellipse with major and minor axes a and b refers to 
the lensing potential rather than the projected density. 
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Fig. 2. — The expected magnification for giant arcs, /it, scaled 
by (rt/W) 2 , as a function of the dimensionless radius s. As in 
Fig. [l] the solid, dotted and dashed lines represent the SIS, NIS, 
and NFW density profiles, respectively. From bottom to top, the 
dashed lines show the results for a NIS profile with r B /r c = 50, 10, 5, 
respectively (see the text for more detail). 

Note that the definition is chosen so as to conserve the 
total projected mass within r G independent of ellipticity, 
as pointed out by Meneghetti et al (2005b). 

Projected ellipticities of real clusters are not very large; 
they are always smaller than their real ellipticities in 
three dimensions. We restrict our theoretical results to 
the first order in e in this subsection. Then, the critical 
curve approximates an ellipse. Without restriction, we 
assume that its major and minor axes are oriented along 
the x and y coordinates, respectively. Giant arcs prefer- 
entially appear near positions with either x « or y w 
from sources in the neighbourhood of cusp points in the 
caustics. 

The gradient of the Jacobian determinant 
(ddet*4/ds)| t of the lens mapping Eq. ((4J perpen- 
dicular to the critical curve at x = and y = 0, where 
giant tangential arcs are expected to appear, is shown 
in Fig. [3] as a function of the dimensionless radius s for 
different ellipticities e. The line types are chosen as in 
Figs. [1] and [5] For a given e, the gradient exhibits the 
same behaviour as in Fig. [T] for axially-symmetric SIS, 
NIS, and NFW density profiles for the same reasons. 

Comparing with axially-symmetric lens models, the 
situation for different ellipticities e is more complicated. 
Two factors will affect the gradient of the Jacobian deter- 
minant (ddet^4/ds)|t and its influence on a giant arc at 
a given position: (1) the total mass enclosed by a circle 
traced by the arc, and (2) the density gradient. Elliptic- 
ity enhances the tidal field and thus shifts critical curves 
towards larger radii. Thereby, it increases the gradient 
(ddet A/ds)\% because of the increasing steepness of the 
projected density profile for a given mass model. 

As in the previous subsection, we plot in Fig. [4] the 
expected magnification for arcs near the x and y axes as 
a function of their distance from the cluster centre for 
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Fig. 3. — Same as Fig.[T] but choosing different ellipticities e ^ 0. 
Giant arcs appearing near the x and y axes are considered (see text 
for further detail). 

asymmetric lens models. The figure shows that the arc 
magnifications show the same general trends as in Fig. [21 
for the same reasons. Note that the magnifications of 
axially-symmetric and axially-asymmetric models differ 
significantly even if the ellipticity e is small. The pre- 
dicted magnifications of arcs for individual mass models 
decrease for increasing e because of similar reasons as 
discussed above, especially for arcs near x = 0, i.e. near 
the major axis. 

In summary, the expected magnifications for giant arcs 
in lenses with different mass models are different even 
though the arc positions can be reproduced by all of 
them. Asymmetry of the mass distribution plays an 
important role for strong lensing. The resulting mag- 
nifications can change substantially even if the system's 
ellipticity is small. Although the absolute magnifications 
for individual arcs or arclets cannot be obtained obser- 
vationally, lens models will be additionally constrained 
by the relative magnifications of arc or arclet pairs from 
the same sources within a lens system. We shall use the 
well-studied clusters A 370 and MS 2137 as illustrations 
in the next section. 

3. ILLUSTRATION: A 370 AND MS 2137 

As discussed in the previous section, different mass 
models are possible for individual galaxy clusters if they 
are only based on the observed positions of giant arcs. 
Relative magnifications of arc or arclet pairs from indi- 
vidual sources in a lens system can be used to constrain 
the mass model further. We will investigate the arc mag- 
nifications in NIS and NFW models, which can well re- 
produce the arc configurations in A 370 and MS 2137, 
respectively. We choose these two clusters for the follow- 
ing reasons: 

• They are prototypical lensing clusters, with A 370 
consisting of two mass components and MS 2137 
having one dominant component. 
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Fig. 4. — Similar to Fig. [2] but adopting different ellipticities e of 
the lensing potential. Results are shown for giant arcs appearing 
near the x and y axes, respectively (see the text for more detail). 

• Both of them have large tangential and radial arcs, 
which mass models are required to fit. 

• Both of them have at least two arc or arclet 
pairs with photometric data available in several 
bands, which allow us to use relative magnifications 
(i.e. flux ratios) of image systems from individual 
sources to further constrain the mass models. 

• High-quality HST images are available for both of 
them, which can be used for detailed strong-lensing 
studies. 

• Both clusters have been thoroughly studied earlier, 
which allows detailed comparisons. 

• For A 370 in particular, we have high-quality im- 
ages in the B, V, R, I, z, J, H and K' bands taken 
by Subaru with a median seeing of ~ 0".5 (Cowie 
ct al. 2008), and in four infrared channels observed 
by IRAC on-board the Spitzer space telescope in 
addition to its HST archive images. The deep op- 
tical and NIR images permit the identification of 
individual images of lensed galaxies in the A370 
region studied in the present paper. The IRAC im- 
ages have > 2" resolution and will be used to study 
in detail the spectral energy distributions (SEDs) 
for the sources combined with the optical and NIR 
data in a forthcoming paper. 

We shall very briefly summarize earlier work on these 
two clusters in the following subsections. Different mass 
models will be applied to the two clusters in order to 
reproduce the main observed features of arcs and arclets. 
Since there are many arclets observed in the field of A 370 
and extensive observational data are available, we will 
study this cluster in more detail. 

Our models improve the previous lens models for 
MS 2137 and A 370 in the following way. Gavazzi et 
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al (2003) studied the density profile for MS 2137 in de- 
tail and Comerford et al (2006) successfully applied the 
NFW density profile to this cluster. Both studies did not 
explicitly consider the relative magnifications between 
image pairs from individual sources as constrained by 
the photometric observations, which we focus on in this 
paper. For A 370, we not only employ more observed 
arclet pairs than Comerford et al (2006) to constrain the 
mass model, but also compare their magnifications with 
the observed photometry. Moreover, we study additional 
arclets in more detail, including their estimated lens red- 
shifts. 

We use the same code as in Comerford et al. (2006) 
to fit mass models to observations. Although only NFW 
profiles were fit to observations there, isothermal density 
profiles are also provided by this code. 

Briefly, it proceeds as follows: Varying the density- 
profile parameters, it strives (1) to make arc sources as 
small as possible, (2) to make arcs and arclets produced 
from these sources match the observed data as closely as 
possible, and (3) to avoid additional images arising from 
the sources or make them faint enough to be compatible 
with the observations. 

These three requirements are quantified by three con- 
tributions to a figur-of-merit function, which is jointly 
minimized to find the best-fitting model parameters. 
This function is not precisely a \ 2 function, as it can- 
not be normalised to the degrees of freedom because the 
intrinsic source size is unknown. See Comerford et al. 
(2006) for more detail. Relative magnifications, i.e. the 
flux ratios, between images in image systems of individ- 
ual sources, are not included in the figure-of-merit func- 
tion when fitting mass models to MS 2137 and A 370 
in the following subsections. They are used in the fur- 
ther discussions on how well the best-fitting models can 
match the photometric observations. 

3.1. MS 2137 

MS 2137 is a rich cD cluster at z = 0.313. Early 
studies by Mellier, Fort & Kneib (1993) and Miralda- 
Escude (1995) found that elliptical isothermal models 
can reproduce the observed arc system remarkably well. 
Miralda-Escude (2002) suggested that its mass distribu- 
tion should be highly elliptic. Comerford et al (2006) 
showed that the cluster can also be modelled with NFW 
components. Based on archival data from the HST and 
the ESO-VLT, Gavazzi et al (2003) and Gavazzi (2005) 
investigated its mass profile in detail and found that it 
can be fit by the isothermal and the NFW density profiles 
alike. Based on a possible fifth image in the very cluster 
centre and on stellar kincmatical data, they concluded 
that the isothermal model may be preferred. 

We fit both the NIS and the NFW mass density pro- 
files to MS 2137, attempting to reproduce the observed 
five arcs including the radial arc. We adopt a single 
lens component at the position of the cD galaxy, simi- 
lar to Comerford et al. (2006). The predicted critical 
curves and caustics of the different best-fitting models 
are shown in Fig. [5l They are almost identical for the 
two types of model, which confirms the results of ear- 
lier studies, since a one-component lens system cannot 
be tightly constrained with arc systems at identical or 
similar redshifts. 

The reproduction of the arc systems by the best-fitting 
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Fig. 5. — Critical curves and caustics for the two best-fitting 
mass models of MS 2137. 
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Fig. 6. — The HST image (upper panel) and the reproduction of 
the observed data (lower panel) by the two best-fitting mass models 
for MS 2137, with (AO, A2, A4) and (Al, A5+A6) being identified 
as counter-images of two different sources. Observed arcs, repro- 
duced images and the reconstructed sources are marked by o, X , 
and *, respectively. The critical curves and caustics of the models 
are marked in black. 

mass models are shown in Fig. [6] adopting the notation 
of Gavazzi et al (2003), where (AO, A2, A4) and (Al, 
A5+A6) are identified as counter-images of two different 
sources. The best-fitting model parameters are given in 
Tab.[T] The figures-of- merit obtained with the complete 
set of input data points read off the observational data 
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TABLE 1 

Parameters of the best-fitting models for MS 2137 and 
A 370, with Column (1): name of the cluster; (2): label of 

LENSES WHERE APPLICABLE; (3): CONVERGENCE AMPLITUDE KS; 
(4) SCALE RADIUS FOR NFW, CORE RADIUS FOR NIS MODELS; (5): 

ELLIPTICITY; (6): POSITION ANGLE; (7): AVERAGED 
FIGURE-OF-MERIT FOR THE TOTAL NUMBER. OF DATA POINTS; (8): 
ADOPTED MASS MODEL. 



cluster 


lens 




r c/s 

(kpc/h) 


e 


e 

(deg) 


x 2 


model 


MS 2137 




0.66 


64 


0.11 


95.7 


0.19 


NFW 






3.12 
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0.11 


95.7 


0.82 


NIS 


A 370 


Gl 


0.16 


254 


0.28 


78.0 


0.16 


NFW 




G2 


0.17 


212 


0.07 


-13.4 




NFW 




Gl 


50.0 


1 


0.20 


75.0 


2.5 


NIS 




G2 


30.0 


1 


0.05 


-10.0 




NIS 



are also listed in the table. It appears that the NFW pro- 
file is preferred based on the figures-of-merit of different 
best-fitting mass models. Note that the figure-of-merit 
in this paper adds the three contributions of sources, the 
observed and the predicted images. Although it does not 
have the usual absolute statistical meaning of a properly 
normalised \ 2 function, it can still serve to evaluate dif- 
ferent models based on identical input data (see Comer- 
ford et al (2006) for more detail). 

Since we use the same multiple- image systems as Com- 
erford et al (2006) to constrain our models, the model 
parameters obtained for the NFW profile are similar ex- 
cept for a very small difference in k s , where we find 0.66 
rather than the 0.67 of Comerford et al (2006), which is 
perfectly within the expected uncertainty. 

The parameter n s is the amplitude of the convergence 
profile k(x), defined by 



k(x) 



2k s 



1 



1 - 



. arctan 



1 



for the NFW profile with x = r/r s and by 

x 2 + 2x1 



k{x) 



2 {x 2 + x 2 ) 3 / 2 



(11) 



(12) 



for the NIS profile, where x and the radius r and 

the core radius r c in units of the Einstein radius. Both 
the NFW and the NIS models predict a fifth-image in the 
central cluster region which was confirmed by Gavazzi et 
al (2003) after subtracting the light produced by the cD 
galaxy. 

The predicted absolute magnifications for the individ- 
ual arcs are shown in Fig. [7] for the two different best- 
fitting mass models for MS 2137, with the solid diagonal 
indicating where the two magnifications for individual 
arcs would agree. As the figure shows, the arc magni- 
fications predicted by the two models agree well. This 
implies that further constraints on the density profile 
for MS 2137 based on relative image magnifications of 
arcs from individual sources cannot be expected because 
MS 2137 is a relatively simple lens. 

Observationally, we can obtain the relative magnifica- 
tions for three arc or arclet pairs, i.e., A0/A2 and A0/A4 
for the giant arc and A1/(A5+A6) for the radial arc, re- 
sulting from two sources based on the photometric data 
in seven bands (U, B, V, R, I, J and K; Gavazzi et 
al 2003). The averaged results obtained in seven bands 
are shown in Fig. [8] with their photometric errors. We 
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Fig. 7. — Predicted magnifications of arcs and arclets in the 
fields of MS 2137 based on the NIS model, compared to those based 
on the NFW model, with (AO, A2, A4) and (Al, A5+A6) being 
identified as image systems of two different sources. The solid line 
illustrates where they would agree. 

find that the model predictions by both the NIS and the 
NFW model are roughly in agreement with the averaged 
relative magnifications for these three pairs within the 
observational errors, although the NFW model again ap- 
pears slightly preferred. We conclude that both the NIS 
and the NFW mass profiles work well for MS 2137, even 
if the magnification, which is the second order of the 
derivative of the potential, is considered. 

Moreover, with the exception of the radial arc, the pre- 
dicted magnifications are somewhat higher for the NIS 
than for the NFW model, which can be understood as 
follows. According to the parameters of the best-fitting 
models for MS 2137 given in Tab. [TJ r s /r c ~ 4, showing 
that the core radius of the NIS model is comparable to 
the scale radius of the NFW profile. This yields higher 
predicted magnifications for the arcs reproduced by the 
NIS compared to the NFW model (cf. Fig.g}. It should 
be pointed out that the predicted magnification of the 
radial arc is calculated according to Eq. {1} based on its 
observed image location, which straddles the correspond- 
ing critical curve. The resulting magnification mainly de- 
pends on its observed size, which is difficult to measure 
accurately since it is located in the very central region of 
the lens. 

3.2. A 370 

Abell 370 (z — 0.375) was the first cluster in which 
gravitational arcs were observed (Soucail et al 1987). 
This two-component, strongly lensing system provides 
an interesting case for us to study the magnifications 
of individual observed arcs and arclets in detail, and to 
compare with previous work. 

This cluster has frequently been studied. For instance, 
Kncib et al (1993) reconstructed the mass distribution 
of A 370 based on the observed giant arc (AO) and the 
arc pair (B2,B3) to conclude that the dark matter ap- 
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Fig. 8. — Comparison of the predicted relative magnifications of 
arcs and arclets pairs in the fields of MS 2137 with the observations. 
Points from left to right denote R(A1)/A5+A6, A0/A2 and A0/A4 
with the photometric errors attached. The solid diagonal indicates 
where they would agree. 

pears to follow the light distribution for this cluster. 
Bartelmann (1996) showed that, as in MS 2137, the ra- 
dial arc can be reproduced by an NFW model, which 
was recently confirmed with a detailed model by Com- 
erford et al (2006). Using a deep HST image, Bezecourt 
et al (1999) constructed a mass model for A 370 based 
on a truncated pseudo-isothermal distribution and found 
that the observed lensed images can be well reproduced. 
Using A 370 as a natural telescope, Hu et al (2002) 
found a galaxy at z w 6.56 and Chary, Stern & Eisen- 
hardt (2005) argued that its star formation rate could 
be > 14OM0yr _1 , thus providing useful clues to the 
galaxy- formation models at high redshift. 

We adopt the notation of Bezecourt et al (1999) and 
Kneib et al (1993) for the different arcs and arclets in 
the field of A 370 to facilitate comparisons with ear- 
lier studies. The additional arclet candidates we identify 
are jointly labeled by "Z" below. For clarity, we only 
show one colour image of A 370 combined from the B, R 
and K' bands rather than separate images in individual 
bands, with the individual arcs and arclets identified in 
Fig. [9] Note that A5 appears near the bottom since it 
appeared near the boundary of the if -band image, but 
it can be very clearly seen in the B band. 

Since there are many arcs and arclets in the field of 
A 370, we list below the image systems produced by in- 
dividual sources for clarity, except the well defined giant 
arc AO, the radial arc R and the arc systems B, C, D and 
E (Bezecourt et al 1999). In the labeled "A" sequence, 
(Al, A2) and (A3, A4) are arc pairs produced by two 
different sources. A5 is a single arc. For the newly found 
arcs labeled by "Z" , we suggest in the next subsection 
that Zl, Z3, and Z6 be arcs produced by three different 
sources and (Z2, Z4, Z5) result from a single source. Z7 
and Z8 (also labeled as W3 and Wl respectively in Fig. 
[9]) with W2 are identified as counter-images of another 




Fig. 9. — Composite image of A 370 taken in the B, R and K' 
bands with the arcs and arclets marked. 



source (see next subsection). 

As for MS 2137, we first construct both a NIS and a 
NFW mass model for A 370 reproducing the giant arc 
AO, the radial arc R, the arclets (Al, A2), (B2, B3), 
and A5. An extra arclet pair labeled (A3, A4) is also 
included, which is easily seen in Fig. [3] especially in the 
-B-band. In total, we use nine relatively bright arcs and 
arclets produced by six sources to constrain the mass 
model for A 370. The lens components are fixed at the 
centres-of-light of the two bright cD galaxies. 

The best-fitting model parameters are listed in Tab. [T] 
together with the averaged figures-of-merit of the com- 
plete set of data points. The corresponding critical 
curves, caustics and the reproduced images are shown 
in Figs. [T0l and ITT1 Our results for the NFW model are 
very close to those found by Comerford et al (2006) where 
the six arcs and arclets AO, R, (Al, A2) and (B2, B3) 
were taken into account. This suggests that (Al, A2) 
and (B2, B3) play important roles in the mass configu- 
ration for A370 (Kneib et al 1993). We point out that 
the high value of k s and the low value of r c in the best- 
fitting NIS model imply that the arc in A 370 can be 
very well reproduced by nearly-singular isothermal lens 
components. 

Since the redshifts of individual sources can also be 
considered as free parameters during the model fitting, 
we can obtain gravitational-lcnsing estimates for indi- 
vidual source redshifts, using the measured redshift of 
the giant arc AO as a reference for the relative lensing 
efficiency at different redshifts. This will be discussed 
further below in this subsection. 

Note that the critical curves and caustics of the best- 
fitting models are very different for the two mass models, 
although both of them reproduce the observed positions 
of individual arcs and arclets very well (cf. Fig. [11]). Since 
the selected arcs and arclets can well constrain the mass 
configuration for A 370, this implies once again that mass 
profiles in clusters cannot be clearly distinguished based 
exclusively on the image positions observed in strong 
lensing. 

Although both best-fitting models succeed in repro- 
ducing the main features of the observed arcs and ar- 
clets, there are four extra images predicted that were 
not yet found observationally. The image near the top 
in Fig.[Tl]is the counter-image predicted for B2 and B3, 
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Fig. 10. — The critical curves and caustics for the best-fitting Fig. 12. — Predicted magnifications of arcs and arclets in the 

mass models for A 370. field of A 370 based on the NIS model, compared to those based 

on the NFW model. The solid diagonal indicates where they would 
agree. 
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Fig. 11. — The best-fitting results for A 370, using the same 
notation as in Fig. [6] 

which matches the observed counter- image labeled B4 in 
the figures of the next subsection. The other three are 
the predicted counter-images of the radial arc. These 
deviations arise because (1) the lens model constructed 
for A 370 only has two components centred on the cD 
galaxies without taking perturbations by additional lens 
components centred on less massive cluster galaxies into 
account; (2) the three additional images are too faint to 
be observed which is consistent with the predictions that 
their magnifications are small. 

The predicted magnifications of the arcs and arclets 
selected to constrain the mass model for A 370 are cal- 
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Fig. 13. — Comparison of the predicted relative magnifications 
of arcs and arclets pairs in the field of A 370 with the observations. 
Points from left to right denote A3/A4 , A1/A2 and B2/B3 with 
the observational errors added, respectively. The solid diagonal 
indicates where they would agree. 

culated for both the NIS and the NFW density profiles 
and displayed in Fig. [12] for comparison. The expected 
magnifications of these nine arcs and arclets are always 
larger for the NFW than for the NIS profile, with the 
magnification of the giant arc being larger by a factor of 
~ 2.5. This is easily understood based on the theoretical 
results shown in Fig.[U Since the ratio r s /r c is very large 
for the best-fitting models of A 370, the magnifications 
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by the NFW profile must exceed those by the N1S profile. 

Since there are many arc and arclet pairs in A370 and 
high-quality images are available in eight optical and N1R 
bands, we can compare the expectations for the relative 
magnifications of the confirmed arclet pairs (B2, B3), 
(Al, A2) and (A3, A4) with the observations in indi- 
vidual bands, to determine which density profile may 
be preferred. Photometric data reduction is done using 
SExtractor (Bertin & Arnouts 1996) for optical and NIR. 
Here, we only show the results, while the detailed pho- 
tometry will be described in a forthcoming paper study- 
ing the SEDs of the lensed galaxies in A 370. 

In Fig. [13j we plot the relative magnifications (i.e., 
flux ratios) of individual arclet pairs averaged over eight 
independent bands against the corresponding predicted 
values. Unlike in MS 2137, we find that the relative mag- 
nifications expected from the NFW profile can match the 
observations much better than those based on the NIS 
profile. The pair (B2, B3) plays an important role here. 
Mellier et al (1991) had already pointed out that (B2, B3) 
is very sensitive to the lens model. We confirm that the 
best-fitting NFW profile achieves a lower figure-of-merit 
than the NIS profile. According to the above analysis 
and current numerical simulations, we thus suggest that 
the NFW profile for A 370 matches the strong-lensing 
observations better than the NIS profile, and adopt it in 
the remainder of this paper. 

As discussed above, the redshifts of individual sources 
can be seen as input parameters in the model fitting be- 
cause lensing by a fixed mass distribution depends on a 
distance factor (see Comerford et al 2006) for sources at 
different redshifts. This allows redshift estimates for in- 
dividual sources required by the best-fitting models. We 
proceed as follows. We choose the giant arc AO as a ref- 
erence, which has the spectroscopic redshift zao = 0.724 
(Soucail et al 1988). By varying the distance factors of 
the other eight arcs and arclets and optimising the model, 
we estimate their redshifts. The error bars are estimated 
by the redshift interval within which the figure-of-merit 
increases by a factor of two from the minimum. 

We find that the source redshifts obtained in this way 
for the best-fitting NFW and NIS models are similar. 
They are shown in the upper part of Table [2] Redshifts 
from spectroscopy, photometry and lensing estimates for 
the same sources obtained in earlier studies are listed 
in Column 4 where available. The table shows that the 
redshifts found here are consistent with the earlier results 
except for the arc pair (Al, A2), which is suggested to 
have a lensing redshift z w 1.4 ± 0.2 (Bezecourt et al 
1999). According to our photometric data in 12 bands, 
we find that Al and A2 are very blue (cf. Fig. [9|) and 
are not detected by IRAC. In qualitative agreement with 
our lensing-based estimate, this suggests (Al, A2) are at 
lower redshift. 

3.3. Further study of the arclets in A 370 

The high-quality and deep images in eight bands from 
B to K' taken by Subaru enable us to apply the best- 
fitting NFW lens model obtained in the preceding subsec- 
tion to some arclets that were already found and labeled 
by Kneib et al (1993) and Bezecourt et al (1999), and 
to some new arclet candidates (labelled "Z") in order 
to test whether the best-fitting model is consistent with 
them (see Fig. [9]). 



TABLE 2 

Redshift estimates z\ (column 3) for arcs in A 370 (column 

1) OBTAINED FROM DISTANCE FACTORS (DF, COLUMN 2) IN THE 

BEST-FITTING LENS MODEL, COMPARED TO SPECTROSCOPIC, 
PHOTOMETRIC OR LENSING- ESTIMATED REDSHIFTS, Z s , Z p OR Z\ 
(COLUMN 4), WHERE AVAILABLE FROM THE LITERATURE (COLUMN 

5). 



arc 


DF 


21 


2 S , z p or zi 


ref. 


AO 


1.00 




0.724 (s) 


[1] 


R 


1.46±0.10 


1.31±0.12 


1.312 (p) 


[2] 








1.7±0.2 (1) 


[3] 


B2,3 


1.10±0.09 


0.80±0.09 


0.806 (s) 


[3] 


Al,2 


1.15±0.11 


0.85±0.11 


1.4±0.2 (1) 


[3] 








> 1-2 (p) 


[-1] 


A3,4 


1.30±0.12 


1.02±0.11 


> 0.2 (p) 


w 


A5 


1.46±0.06 


1.30±0.12 


1.306 (s) 


[5] 


B4 


1.11 


0.80 


- 


- 


Cl,2,4 


1.05±0.13 


0.76±0.12 


0.75±0.1 (1) 


m 








0.60 - 0.92 (p) 


w 


Dl,2 


1.10±0.11 


0.80±0.11 


0.85±0.1 (1) 


[3] 








< 0.4 (p) 


[-1] 


El,2,3 


1.30±0.11 


1.02±0.18 


1.3±Q.l (1) 


[3] 


Zl 


1.70±0.05 


2.28±0.29 






Z2 


0.95±0.14 


0.69±0.10 






Z3 


1.60±0.06 


1.73±0.21 






Z4 


0.95±0.15 


0.69±0.11 






Z5 


0.95±0.14 


0.69±0.10 






Z6 


l.OOiO.ll 


0.72±0.09 






Z7 


1.34±0.08 


1.08±0.11 






Z8 


1.34±0.08 


1.08±0.11 






W2 


1.34±0.08 


1.08±0.11 







DF: distance factor (Comerford et al 2006); z\\ lens redshift; z B : 
spectroscopic redshift; z p : photometric redshift; (s): spectroscopy; 
(p): photometry; (1): lens. 

Refs: [1] Soucail et al (1988); [2] Smail et al (1996); [3] Bezecourt 
et al (1999); [4] Kneib et al (1994); [5] Mellier et al (1991) 

We thus apply the best-fitting NFW lens model ob- 
tained in the previous subsection to the arclet pair candi- 
dates labeled B, C, D and E by first optimising the model 
under the assumption of minimal source sizes, only us- 
ing their redshifts as free parameters. The results of this 
fit for individual arclet pairs are shown in Fig. 1141 and 
their estimated redshifts are listed in the middle part of 
Tabled 

Figure [T4l shows that the best-fitting NFW model can 
fairly reproduce the observed images of the arclet pairs 
B, C, and D. Note that the additional image produced by 
(B2, B3) in Fig. [11] matches well their possible counter 
image B4. Their estimated redshifts are consistent with 
previous photometric measurements. Small differences of 
the image positions between the model predictions and 
the observations for E3 in the right panel of the figure 
may arise because we ignore perturbations due to lensing 
by less luminous galaxies within A 370. Such differences 
in the image positions measured and reproduced by the 
model also occur for C4. For the pair (Dl, D2), the 
model predicts an extra counter-image very close to Dl 
that is not clearly observed, possibly because the image 
was yet too faint to be seen. 

In addition to the arclet candidates labeled by Kneib 
et al (1993) and Bezecourt et al (1999), we identify new 
arclet candidates in the field of A 370 labeled "Zl" to 
"Z8" and shown in Fig. [9) with the measured ellipticities 
of individual objects larger than 0.6 in R band by SEx- 
tractor (Bertin & Arnouts 1996) and checking by eye in 
the other bands. Similar to B, C, D and E, the best- 
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Fig. 14. — The best-fitting results for the arclet pairs B, C, D, 
E in A 370. 
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in A 370. 



The best-fitting results for the new arclet candidates 



fitting NFW model is also applied to them, taking their 
rcdshifts as free parameters. Best-fitting results for them 
are shown in Fig. [15] with the expected lensing redshifts 
for individual new arclets being listed in the final part of 
Table[2j A detailed study of their SEDs will be presented 
in a separate paper. 

Based on these results for the new arclets, we suggest 
that Z6 may be a lensed galaxy at low redshift. Z2, Z4 
and Z5 could be counter-images because their recovered 
sources are close to each other in the sky, and they have 
similar lensing redshifts. Although the predicted sepa- 
rations of the source positions between Z4 and Z2, and 



between Z4 and Z5 are possibly too large, this may be 
due to model uncertainties. Moreover, two candidates of 
distant objects, Zl and Z3, are found with lensing red- 
shifts of ~ 2.3 and ~ 1.73, respectively, which need to be 
spectroscopically or photometrically confirmed. 

Finally, the new candidates Z7 and Z8 are suggested to 
be counter-images of a source at redshift near 1.3. Based 
on the best-fitting NFW model, we also found their pre- 
dicted counter image close to B2, labeled W2 in Figs. [TBI 
and [5] Photometry confirms that it has the same colour 
as Z7 and Z8. We note that this newly found arclet 
pair had already been detected by Bernard Fort and col- 
laborators several years ago, but remained unpublished 
(B. Fort, private communication). Like (B2, B3), these 
counter images are also very sensitive to the mass model 
of A 370. Since they are extremely red, the study of their 
SEDs in the next paper suggests that they arise from an 
active star- forming galaxy heavily enshrouded by dust. 

4. CONCLUSIONS 

Since it is virtually impossible to distinguish mass 
models with different density profiles for galaxy clusters 
based exclusively on the observed arc and arclet posi- 
tions, which depend on the first derivative of the effec- 
tive lensing potential, we have studied in detail the im- 
age magnifications expected from lens models assembled 
from components with isothermal and NFW density pro- 
files. We considered both axially-symmetric and asym- 
metric models and conclude that the magnifications for 
different mass models can be substantially different. The 
ellipticity of a mass model must be taken into account 
because it plays an indispensable role in strong-lcnsing 
studies. 

Given a lens model reproducing the position of a gi- 
ant arc, its predicted magnification is usually smaller 
for singular than for non-singular isothermal profiles. 
Isothermal models with small cores predict magnifica- 
tions smaller than models with NFW mass profiles. 
Isothermal mass models with large cores, whose density 
gradient in the inner region is nearly flat, predict mag- 
nifications increasing with the core radius and exceeding 
those for NFW density profiles. These effects are more 
pronounced for increasingly asymmetric models. 

We illustrate our results with the two well-studied clus- 
ters MS 2137 and A 370, which we model using the code 
described by Comerford et al (2006). Comparing with 
previous studies, the main improvements in this paper 
are based on (1) the photometry for obtaining relative 
magnifications of individual counter-image pairs to con- 
strain the lens models; and (2) further studies of addi- 
tional arclets in A 370 including lens redshifts of newly 
found arclets. 

The observational data were taken with Subaru, HST 
and Spitzer for A 370 and taken from Gavazzi et al 
(2003) for MS 2137. We find that both isothermal and 
NFW models can reproduce the observed arcs and ar- 
clet positions in MS 2137 and A 370, which is consis- 
tent with the results of Comerford et al (2006), Gavazzi 
(2005), Gavazzi et al (2003), Miralda-Escude (2002, 
1995), Bezecourt et al (1999), Mellier, Fort & Kneib 
(1993), and Kneib et al (1993). 

Relative magnifications of individual arc and arclet 
pairs are adopted to adjust the mass model by compar- 
ing with observations. Since MS 2137 is a relatively sim- 



Comparisons of mass profiles in strong lensing 



11 



pie, single-component system, it is hard to distinguish 
isothermal from NFW mass models. On the other hand, 
A 370 has two major lens components with several arc 
and arclet pairs, which allow us to adopt the relative 
magnifications of the confirmed arclet pairs to conclude 
that the NFW profile is more realistic than the isother- 
mal profile for A 370. 

The best-fitting NFW mass model for A 370 is applied 
to several arclet candidates found by Kneib et al (1993) 
and Bezecourt et al (1999) that were not used for con- 
straining the mass model. We find that this best-fitting 
model can well reproduce the observed arclets. 

The high-quality images of A 370 in twelve bands from 
B to infrared will allow us to recover SEDs for individ- 
ual sources in the field of A 370 in a forthcoming paper. 
We also find eight new arclet candidates, apply our best- 
fitting NFW lens model to them and show that we can 
reproduce their observed features fairly well. We ob- 
tain redshift estimates from lensing for individual arcs 
and arclets in A 370 that are consistent with spectro- 
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scopic and photometric results where available from ear- 
lier studies. Other redshift predictions from lensing need 
to be confirmed by future observations. Finally, we found 
a new arclet pair labeled Z7, Z8 and W2 and suggest that 
they arise from an active star-forming galaxy heavily en- 
shrouded by dust. 
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